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Abstract. We study the ferromagnetic phase transition in a randomly layered Heisenberg 

T-r^ , model. A recent strong-disorder renormalization group approach [Phys. Rev. B 81, 144407 

r^ • (2010)] predicted that the critical point in this system is of exotic infinite-randomness type and is 

Q ' accompanied by strong power-law Griffiths singularities. Here, we report results of Monte-Carlo 

O . simulations that provide numerical evidence in support of these predictions. Specifically, we 

investigate the finite-size scaling behavior of the magnetic susceptibility which is characterized 
by a non-universal power-law divergence in the Griffiths phase. In addition, we calculate the 
►^ ■ time autocorrelation function of the spins. It features a very slow decay in the Griffiths phase, 

^^ ' following a non-universal power law in time. 

in 
in 

f^ ■ 1. Introduction 

^D . Impurities, defects, and other types of quenched disorder influence zero-temperature quantum 

phase transitions much more strongly than generic classical phase transitions that occur at 
nonzero temperatures. The interplay of quantum and disorder fluctuations produces uncon- 

k> ', ventional phenomena such as power-law quantum Griffiths singularities [1, 2, 3], infinite- 

$H ' randomness critical points featuring exponential instead of power-law scaling [4, 5], and smeared 

phase transitions [6, 7]. A recent review of these phenomena can be found in Ref. [8], while Ref. 
[9] focuses on metallic systems and also discusses experiments. 

Disorder effects are stronger at quantum phase transitions than at classical transitions because 
the quenched disorder is perfectly correlated in the imaginary time direction. Imaginary time 
acts as an extra dimension at a quantum phase transition and becomes infinitely extended at 
zero temperature. Therefore, the impurities and defects are effectively "infinitely large" in this 
extra dimension, which makes them much harder to average out than finite-size defects. 

This argument suggests that strong disorder effects should also occur at a classical thermal 
phase transition provided that the disorder is perfectly correlated in one or more space 
dimensions. An example occurs in the McCoy- Wu model, a disordered classical two-dimensional 
Ising model having perfect disorder correlations in one of the two dimensions. McCoy and Wu 
[10, 11, 12, 13] showed that this model exhibits an unusual phase transition featuring a smooth 
specific heat while the susceptibility is infinite over an entire temperature range. Fisher [4, 5] 
achieved an essentially complete understanding of this phase transition with the help of a strong- 
disorder renormalization group approach (using the equivalence between the McCoy- Wu model 
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and the random transverse- field Ising chain). He determined that the critical point is of exotic 
infinite-randomness type and is accompanied by power-law Griffiths singularities. 

Recently, some of us investigated another classical system with perfectly correlated disorder, 
the randomly layered Heisenberg magnet, by means of a strong-disorder renormalization group 
[14]. This theory predicts that the (thermal) critical point in the randomly layered Heisenberg 
magnet is of infinite-randomness type as well. Moreover, it is in the same universality class as 
the quantum critical point of the random transverse-field Ising chain. In this paper, we present 
the results of Monte-Carlo simulations of the randomly layered Heisenberg model. They provide 
numerical evidence in support of the above renormalization group predictions. 

2. Model and renormalization group predictions 

We consider a ferromagnet consisting of a random sequence of layers made up of two different 
ferromagnetic materials. Its Hamiltonian, a classical Heisenberg model on a three-dimensional 
lattice of perpendicular size L±_ (in z direction) and in-plane size Lh (in the x and y directions) 
is given by 

H = -2_^Jl (Sr • Sr_|_x -I- Sr • Sp+y) - 2^ J z Sr • Sr+z- (1) 

r r 

Here, Sr is a three-component unit vector on lattice site r, and x, y, and z are the unit vectors 
in the coordinate directions. The interactions within the layers, JI, and between the layers, 
J^, are both positive and independent random functions of the perpendicular coordinate z. In 
the following, we take all J^ to be identical, J^ = J-*-, while the JI are drawn from a binary 
probability distribution P{J^^) = (1 — p) (5(jll — J„) +p5{j\\ — Ji) with J„ > Ji. Here, p is the 
concentration of the "weak" layers. 

The qualitative behavior of (1) is easily explained. At sufficiently high temperatures, the 
model is in a conventional paramagnetic phase. Below a temperature T„ (the transition 
temperature of a hypothetical system having JI = Ju for all z) but above the actual critical 
temperature Tc, rare thick slabs of strong layers develop local order while the bulk system is still 
nonmagnetic. This is the paramagnetic Griffiths phase (or Griffiths region). In the ferromagnetic 
Griffiths phase, located between Tc and a temperature Ti (the transition temperature of a 

hypothetical system having Jz = Ji for all z), bulk magnetism coexists with rare nonmagnetic 
slabs. Finally, below T;, the system is in a conventional ferromagnetic phase. 

In Ref. [14], the behavior in both Griffiths phases and at criticality has been derived within 
a strong-disorder renormalization group calculation. Here, we simply motivate and summarize 
the results. The probability of finding a slab of -L_r_r consecutive strong layers is given by simple 
combinatorics; it reads w{LjiJi) '-^ (1 — p^^RR = ^-p^rr with p = — ln(l — p). Each such 
slab is equivalent to a two-dimensional Heisenberg model with an effective interaction LjinJu- 
Because the two-dimensional Heisenberg model is exactly at its lower critical dimension, the 
renormalized distance from criticality, e, of such a slab decreases exponentially with its thickness, 
e{Ljiji) '-^ e~^^^^ [8, 15]. Combining the two exponentials gives a power law spectrum of locally 
ordered slabs, 

p(e) ~ eP/"-^ = €^l''-^ (2) 

where the second equality defines the conventionally used dynamical exponent, z. It increases 
with decreasing temperature throughout the Griffiths phase and diverges as z ^ 1/|T — Td at 
the actual critical point. 

Many important observables follow from appropriate integrals of (2). The susceptibility can 
be estimated by x ~ I de p{e)/e- In an infinite system, the lower bound of the integral is 0; 
therefore, the susceptibility diverges in the entire temperature region where z > 1. A finite 
system size L\\ in the in-plane directions introduces a nonzero lower bound emm ~ L7 . Thus, 
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Figure 1. Left: Susceptibility x ^-s a function of in-plane system size Lm for several temperatures 
in the Griffiths region. The perpendicular size is L± = 800; the data are averages over 160 
disorder configurations. The solid lines are fits to the power law (3). Right: Resulting Griffiths 
dynamical exponent z as a function of temperature T in the weakly disordered Griffiths region. 



the susceptibility in the Griffiths region diverges as 
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for z > 1. The behavior of the time autocorrelation function C{t) = (1/A^) X]r('S(r, i) • S(r, 0)) 
within model A dynamics [16] can be determined analogously {N is the total number of 
spins). As the correlation time of a single, locally ordered slab is inversely proportional to 
its renormalized distance from criticality, we obtain, in the Griffiths phase. 
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3. Monte-Carlo simulations 

We have carried out Monte-Carlo simulations of the Hamiltonian (1) with system sizes up to 
L_|_ = 800 and Lu = 150. While studying the thermodynamics, we have used the efficient Wolff 
cluster algorithm [17] to eliminate critical slowing down. To measure the time autocorrelation 
function, we have equilibrated the system using the Wolff algorithm but then propagated the 
dynamics by means of the Metropolis algorithm [18] which implements model A dynamics. All 
results reported below are for J^ = 1.00 and a concentration p = 0.8 of weak layers with 
Ju = 1.00 and J; = 0.25. The data are averages over a large number (20 to 160) of disorder 
realizations. 

To test the finite-size behavior (3) of the susceptibility, one needs to consider samples having 
sizes L± ^ Lii such that L± is effectively infinite. We have used system sizes L± = 800 
and Lii = 35 to 100. Figure 1 shows the susceptibility x as a function of Lu for several 
temperatures in the Griffiths region below T„ w 1.443. x follows a nonuniversal power law 
in Lii with a temperature-dependent exponent. Simulations for many more temperature values 
yield analogous results. The value of the exponent z extracted from fits to (3) is shown in the 
right panel of figure 1 for the paramagnetic side of the Griffiths region, z can be fitted to the 
predicted power law z ~ 1/|L — Tc\, as discussed after (2), giving the estimate Tc w 0.933 

Figure 2 shows the autocorrelation function C{t) of the spins as a function of time t for 
several temperatures in the Griffiths phase. In agreement with the prediction (4), the long- 
time behavior of C{t) follows a nonuniversal power law. We have also studied (not shown) the 
behavior at criticality. It displays the logarithmic time-dependence predicted in Ref. [14]. 




Figure 2. Time autocorrelation function 
C(t) for several temperatures in the 
Griffiths phase. The system sizes are 
L_L = 400 and Ly = 100. The data are 
averages over 80 disorder configurations. 
The solid lines are fits to (4). 



4. Conclusions 

In summary, we have performed Monte-Carlo simulations of a randomly layered three- 
dimensional Heisenberg model. Our numerical results for the magnetic susceptibility and 
the time autocorrelation function in the paramagnetic Griffiths phase support the infinite- 
randomness critical point scenario that arises from the strong-disorder renormalization group 
approach [14]. A full confirmation of the theoretical predictions requires performing a scaling 
analysis of the infinite-randomness critical point itself, and measuring all three independent 
critical exponents. This work is in progress. 

Experimental verifications of infinite-randomness critical behavior and the accompanying 
power-law Griffiths singularities have been hard to come by, in particular in higher-dimensional 
systems. Only very recently, promising measurements have been reported [19, 20] of the quantum 
phase transitions in CePdi-^iRh^; and Nii-^Va;. The randomly layered Heisenberg magnet 
considered here provides an alternative realization of an infinite-randomness critical point. It 
may be more easily realizable in experiment because the critical point is classical, and samples 
can be produced by depositing random layers of two different ferromagnetic materials. 

This work has been supported in part by the NSF under grant nos. DMR-0339147 and 
DMR-0906566 and by Research Corporation. 
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